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We study a spinless level that hybridizes with a fermionic band and is also coupled via its charge 
to a dissipative bosonic bath. We consider the general case of a power-law hybridization function 
r(w) oc \ui\ r with r > 0, and a bosonic bath spectral function B(ui) oc cu s with s > —1. For 
r < 1 and max(0, 2r — 1) < s < 1, this Bose-Fermi quantum impurity model features a continuous 
zero-temperature transition between a delocalized phase, with tunneling between the impurity level 
and the band, and a localized phase, in which dissipation suppresses tunneling in the low-energy 
limit. The phase diagram and the critical behavior of the model are elucidated using perturbative 
and numerical renormalization-group techniques, between which there is excellent agreement in the 
appropriate regimes. For r = this model's critical properties coincide with those of the spin-boson 
and Ising Bose-Fermi Kondo models, as expected from bosonization. 



I. INTRODUCTION 

Quantum phase transitions^ in mesoscopic systems 
form a growing area of condensed matter research. From 
a theoretical perspective, it is known that models of 
a finite system (the "impurity") coupled to infinite 
baths may exhibit boundary quantum phase transitions 
(QPTs), at which only a subset of the degrees of free- 
dom becomes critical^ Such models help to advance our 
understanding of quantum criticality in strongly corre- 
lated systems: Concepts and solution techniques devel- 
oped in the impurity context may be applied to lattice 
models, e.g., within the framework of dynamical mean- 
field theory (DMFT)£ and its extensions. This approach 
has been followed in connection with the "local critical- 
ity" proposed to underlie the anomalous non- Fermi-liquid 
behavior of several heavy-fermion systems^ On the ex- 
perimental side, QPTs in mesoscopic few-level systems 
are of great interest, both for the unprecedented oppor- 
tunity to probe quantum criticality in a direct and highly 
controlled fashion^ and for their numerous potential 
technological applications, e.g., in nanoclcctronics and 
quantum information processingf 7 -^ 

In recent years, QPTs have been identified and studied 
in a number of quantum impurity models. 2 Such mod- 
els can contain both fermionic bands (e.g., conduction- 
electron quasiparticles) and bosonic baths (e.g., phonons, 
spin fluctuations, or electromagnetic noise). Analytical 
and numerical techniques have been refined to analyze 
the critical behavior of these models. Analytical ap- 
proaches based on bosonization or conformal field theory 
have been used extensively, although their applicability is 
limited, e.g., to certain forms of the bath spectrum. For 
other situations, powerful epsilon-expansion techniques 
have been developed. As such expansions are asymp- 
totic in character, a comparison with numerical results is 



mandatory to assess their reliability. 

An example with especially rich behavior is the 
fermionic pseudogap Kondo model, 10 which features 
QPTs between Kondo-screened and local-moment 
ground states. -^i 11 ' 12 ' 13 ' 14 ' 15 Essentially perfect agree- 
ment between the results of various epsilon expan- 
sions (around different critical dimensions) and numer- 
ical renormalization-group (NRG) calculations has been 
found in critical exponents as well as universal ampli- 
tudes such as the residual impurity entropy. 14 i 15 

Impurity models that include bosons are harder to 
tackle numerically than pure-fermionic problems due to 
the large Hilbert space, and fewer results are available. 
The development of a bosonic versio n 16 ' 17 of Wilson's 
NRG approach! 8 - has made possible a detailed nonper- 
turbative study of the spin-boson model, where tunnel- 
ing in a two-state system competes with dissipation. 19 
For the case of Ohmic dissipation, the spin-boson model 
has long been known to display a QPT of the Kosterlitz- 
Thouless type. In the sub-Ohmic case, the model instead 
exhibits a line of continuous QPTs governed by interact- 
ing quantum critical points (QCPs). 16 ' 17 ! 20 (The latter 
lie in a different universality class than the QCP of the 
pseudogap Kondo model.) 

Of particular interest, both for mesoscopics and in 
the context of extended DMFT for correlated lattice- 
systems (SIi 2 ^ are impurity models with fermionic and 
bosonic baths. The best-studied member of this class 
is the Bose-Fermi Kondo model ) 23 ' 24 ' 25 ' 26 ' 27 with a spin- 
i local moment coupled to fermionic quasiparticles (the 
regular Kondo model) as well as to a bosonic bath. The 
latter may describe spin or charge fluctuations of the bulk 
system in which the impurity is embedded. The scope of 
NRG applications has recently been widened to provide 
a comprehensive treatment of an Ising-symmetric version 
of the Bose-Fermi Kondo model. 28 ' 29 

The purpose of this paper is to investigate a some- 
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what simpler quantum impurity model containing both 
fcrmionic and bosonic baths, namely a resonant-level 
model of spinless electrons, with the impurity charge cou- 
pled to a dissipative reservoir. In standard notation, its 
Hamiltonian is 



« = e// t / + Z>*(/ t ° fc + H.c.) +J2 



J. 



(/V-DE^+O+E 



u q b q b q . 



(1) 



with Vk characterizing the hybridization between conduc- 
tion electrons of energy Ek and the impurity level at en- 
ergy e f , and g q coupling bosons of energy w q to the im- 
purity occupancy. Without loss of generality, Vk and g q 
are taken to be real and non-negative. Equation {T]) rep- 
resents perhaps the simplest nontrivial Bose-Fermi quan- 
tum impurity model, making it a paradigm for this class 
and an ideal problem for detailed comparison between 
analytical and numerical results. 

The model is completely specified by the impurity level 
energy e f , the hybridization function 



T(w) = ttE «fc<5( w -£fc) = r o 



for 



<D, (2) 



and the bosonic bath spectral function 



B{u) = *Tgl5(u-uj q ) = B Q ( — 



for < u> < uj r 



(3) 



with D and oj c acting as fermionic and bosonic cutoffs, re- 
spectively. Thus, in addition to a power-law spectrum for 
the bosonic bath density of states (DOS) characterized 
by an exponent s, we consider a nonconstant particle-hole 
(p-h) symmetric hybridization function characterized by 
an exponent r. Increasing r (and hence depleting the hy- 
bridization function around the Fermi level u> = 0) and 
increasing Bo both act to suppress tunneling between 
the local level and the conduction band. For most of the 
numerical work presented in Sec. IIII1 we fix r, s, and 
the hybridization strength To, then tune the dissipation 
strength Bq to the vicinity of a QPT. 

Although the bath densities of states and Vk, g q do 
not require separate specification, it will facilitate com- 
parison between numerical and perturbative results to 
assume that Vk — vq, g q — go for all k, q. In this 
case, T(u>) = ttvqP c (uj) and B(oS) — -KffoPbijjj), with the 
fermionic and bosonic DOS given, respectively, by 



p c {u) = iV |«/£>| r 
p b {io) = (K^/tt) (u/u c y 



for|w|<£>, (4) 
for < uj < cu c , (5) 



where No and Kq are normalization factors. Thus, To = 
ttNoVq and B = (K go) 2 - The metallic case is recovered 
for r = 0, and Ohmic dissipation corresponds to taking 
s = 1. 



Delocalized l-2r+s=0 
(Simp- rln 2) 



QPT 



Localized 
(Sim P = In 2) 



FIG. 1: Schematic phase diagram of the dissipative resonant- 
level model (fl| , in the parameter space spanned by exponents 
r and s characterizing the low-energy behavior of fermionic 
and bosonic baths, respectively. (A finite coupling to both 
baths is assumed.) For max(0, 2r — 1) < s < 1, the model 
shows a boundary quantum phase transition (as the couplings 
t>o and go are varied) between a delocalized phase and a lo- 
calized phase. (The physics along the line r = is identical 
to that of the spin-boson model.) In contrast, for s > 1 and 
r < 1 the system is generically delocalized, whereas it is lo- 
calized in the rest of the parameter regime (shaded). The 
impurity entropy Simp is discussed in the text. Two pertur- 
bative RG expansions are employed: around the free-impurity 
fixed point, where the expansion is controlled about r = s = 1 
(Sec. Ill B|) and around the resonant-level fixed point, where 
the expansion is controlled in 1 — 2r + s fSec. Ill C)) . 



It is convenient to identify a pseudospin — making clear 
the close relationship between model |T]) and the spin- 
boson model and its variants — by writing 



f = S+, f = s~ 



s z 



(6) 



In the model described by Eq. |T]), the friction caused 
by the bosonic bath competes with the resonant tunnel- 
ing of electrons. In contrast to the simpler spin-boson 
modelj^ the tunneling properties are determined by the 
hybridization function T(co). 

For Ef — the model features a Z2 symmetry of 
particle- hole type [assuming p e (w) — p c {—Lo) as noted 



above], namely Ck 



k- 



f 



-P, and S z 



Then, we expect that the competition between resonant 
tunneling and dissipation yields a QPT between a "de- 
localized" phase ((S z ) — 0), in which the principal effect 
of dissipation is to renormalize the tunneling amplitude, 
and a "localized" phase ((S z ) ^ 0) with a doubly degen- 
erate ground state, where the tunneling amplitude renor- 
malizes to zero in the low-energy limit. We note that for 
the case of a metallic fermionic bath [r — in Eq. Q], 
bosonization techniques can be used to map the model 
([T|) to the spin-boson model, 3 2 (The same applies to the 
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Ising-symmetric Bose- Fermi Kondo model with r = 0, 
and this equivalence has been verified using NRG. 28 ' 29 ) 
In this paper, we employ renormalization-group (RG) 
techniques to map out the phase diagram of the Hamil- 
tonian (JXJ) and to establish over what range of bath ex- 
ponents r and s the model can be tuned to a delocalized- 
to-localized QPT, akin to that of the spin-boson model. 
We do so using both perturbative RG methods, based 
on epsilon-expansion techniques developed in the con- 
text of the pseudogap Kondo and Anderson models ji 5 - 
and the Bose-Fermi extensio n 28 ' 29 of the NRG approach, 
which allows us to access the entire parameter range of 
the model. 

Our main result is summarized in Fig. [TJ which illus- 
trates the qualitative behavior of the model in the plane 
spanned by the bath exponents r and s. A delocalized-to- 
localizcd transition — which for r = is identically that 
of the spin-boson model — is present at r > as well. A 
more detailed discussion is given in Sec. Ill C 31 

The remainder of the paper is organized as follows. 
The perturbative RG analysis is outlined in Sec. [Til where 
results for various critical exponents are obtained by ex- 
pansion around two distinct fixed points. In Sec. IIII1 we 
provide nonperturbative NRG results for the model, in- 
cluding discussion of the phase diagram, the response to 
a local field, and the single-particle spectral function. We 
find excellent quantitative agreement between analytical 
and numerical results in the appropriate limits. Although 
the critical properties of the model flj) for r = are es- 
tablished via the mapping to the spin-boson model, we 
confirm the equivalence by direct calculation. 



at the FImp fixed point for s < 1 (see, e.g., Refs. [l9l 
and [H). It tends to suppress tunneling in the low-energy 
limit. By analogy with the spin-boson model, this can be 
expected to result in a doubly degenerate ground state, 
Shnp = In 2, i.e., a phase with broken Z2 symmetry. This 
localized fixed point (Loc) corresponds to coupling values 
(g,v) — (oo,0). (Note that for s > 1 the effect of the 
bosonic bath is weak, not causing localization.) 

The preceding discussion suggests that, for r < 1 and 
s < 1, a QPT separates a delocalized (small-dissipation) 
phase from a localized (large-dissipation) phase. Clearly, 
this applies only to the case of p-h symmetry, Sf = 0. 
Otherwise the Z<i symmetry of the Hamiltonian is bro- 
ken from the outset, and the phase transition upon vari- 
ation of the dissipation strength will be smeared into a 
crossover; this is analogous to the behavior of the spin- 
boson model in the presence of a finite bias. Furthermore, 
in situations where the system is localized at et = 0, 
there will be a first-order transition upon tuning e / from 
positive to negative values (as in an ordered magnet sub- 
ject to a field). 

We now proceed with an RG treatment of the model 
(JTJ, carried out without recourse to bosonization. We 
can access quantum-critical properties via two distinct 
expansions: (i) an expansion around the free- impurity 
fixed point (Sec. Ill Bp . which is formally valid provided 
that the couplings to both baths are small, and (ii) an ex- 
pansion around the resonant-level fixed point (Sec. Ill C]) . 
performed after exactly integrating out the c fermions. 
The second approach proves to have the wider range of 
applicability. 



II. PERTURBATIVE RENORMALIZATION 
GROUP 

A. Zero-temperature phases 

We begin by discussing the trivial fixed points of the 
model ([1} in the presence of p-h symmetry, e$ = 0. As a 
characterization, we will refer to the residual impurity en- 
tropy <Si m p, which is defined as the impurity contribution 
to the total entropy in the limit temperature T — > 0. 2 

For vq = go = 0, the impurity is decoupled from both 
baths. We denote this free-impurity fixed point by FImp. 
The ground state is doubly degenerate: Simp = In 2. 

For vo 7^ and go = one has a resonant-level model 
with a power-law conduction-band DOS given by Eq. 
(H]). The hybridization is relevant in the RG sense (w.r.t. 
FImp) for r < 1, and hence the impurity charge strongly 
fluctuates ! 12 ' 15 We refer to this as the delocalized fixed 
point (Deloc), which, as discussed in Ref. HH, is located 
at intermediate RG coupling, (g,v) = (0, «*). Somewhat 
surprisingly, the impurity entropy is <Si mp = rln2, and 
vanishes only in the metallic case r = 0. For r > 1, 
by contrast, the hybridization is RG- irrelevant, and the 
delocalized fixed point merges with Flmp^i 

The dissipative coupling go turns out to be RG-relevant 



B. RG expansion around the free-impurity limit 

In this subsection, we apply an RG epsilon expan- 
sion for weak couplings near the free- impurity fixed point 
where vq = go = 0. 



1. RG equations 

We model the bosonic bath by a relativistic scalar field, 
<j> = b + w , in d = 2 + s dimensions, with the action 



,/3 

S<p = / dr 
Jo 



" d d q 
(2n) d 



0_ <J (r)(-a 2 + (7 2 )^(r), (7) 



A 9 being a momentum-space cutoff (related to the energy 
cutoff co c via to c — cA q with c = 1 being a velocity). This 
produces a DOS of the form 



= sgnH^-^M 8 = sgn(w)^2- 

Z 7T 



(8) 



for \u\ < uj c , with S d = 2/[(47r) d / 2 r(d/2)]. [Note that 
p<j, is just a symmetrized version of pb defined in Eq. 
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(a) 



(b) 



O 



-o 



(c) 
/ 



FIG. 2: Diagrams appearing in the perturbative expansion 
for the dissipative resonant-level model. Dashed, solid, and 
wiggly lines denote respectively /, c, and <j> propagators. The 
gray (black) circles are the interaction vertices v (<?). (a) and 
(b): / fermion self-energy diagrams to one- loop order, (c) 
One-loop vertex renormalization of g. 



©.] Similarly, we represent the fermionic bath by Dirac 
fermions in (1 + r) dimensions: 



S c = I dr 



A,, 



dk\k\ r 



(2< 



l+r 



c k (d T + k)c k , 



(9) 



with Afe = D/vf and vf — 1 being the (Fermi) velocity, 
which reproduces the DOS defined in Eq. (UJ). A path- 
integral representation of Eq. Jl} reads 

S = S C + S^+ drf d T f + go / dr (// - |) cj>(r, 0) 



vo 



/ dr [/c(r, 0) + c.c] . 
Jo 



(10) 



Power counting yields the bare scaling dimensions of 
fields and couplings with respect to vo = go = 0: [/] = 0, 
[<j> q ] = -(1 + s)/2, [c k ] = -(1 + r)/2, [vo] = (1 - r)/2, 
and [go] = (1 — s)/2. Thus, we can carry out an RG 
expansion around r — 1 and s = 1, where both vq and go 
become marginal, defining 



(11) 



2 V " 2 

In order to proceed with the RG analysis, we define a 
renormalized field and couplings v and g according to 



v 



.9o 




(12) 



KoZ f 



where /x is an arbitrary renormalization energy scale and 
Zf, Z v , and Z g are renormalization factors. As is usual 
for impurity problems, there is no renormalization of the 
bosonic and fermionic bulk propagators, since the im- 
purity only provides a one-over-volume correction to the 
bulk properties. The relevant diagrams for obtaining the 
one-loop RG beta functions are shown in Fig. [2] 

Following standard procedures^ the one-loop RG 
beta functions of the dissipative resonant-level model are 
given by 



(3{v) = -e'v + v 3 + \g 2 v, 



(13) 



where the calculation parallels that of Ref. [l5|. The cor- 
responding Z factors, to one-loop accuracy, are Zf = 
1 - v 2 /e' - 5 2 /2e, Z v = 1, and Z g = 1 - g 2 /2e. 

The RG flows arising from Eqs. (fRl)) are plotted in Fig. 
[3] In this subsection, we consider the case < s < 1; the 
regime s < is discussed in Sec. Ill C 21 Fixed points at 



iff 



2 ) = (0,e') and (g* 2 ,v* 2 ) = (oo,0) describe the 



delocalized (Deloc) and localized (Loc) phases, respec- 
tively. For r < r+, where 



r+ = (1 + s)/2, 



(14) 



both these fixed points are stable: For small go and 
large vo, the ground state is delocalized, characterized 
by strong local charge fluctuations due to resonant tun- 
neling between the impurity and the conduction electron 
bath ((S z ) = 0). In the opposite limit of small vo and 
large go, we find a localized ground state where charge 
tunneling renormalizes to zero in the low-energy limit 
((S z ) 7^ 0). An unstable critical fixed point [Cr], located 
at (g* 2 ,v* 2 ) = (2e' - e,e/2), controls the QPT between 
these two phases. This critical fixed point lies on the sep- 
aratrix specifying the phase boundary in the go~vo plane 
between the delocalized and localized phases. 

As r approaches r + from below, the critical fixed 
point merges with the delocalized fixed point (which it- 
self merges with FImp as r — > 1 from below). Hence, no 
transition occurs for r > r + : Deloc and FImp are unsta- 
ble w.r.t. infinitesimal bosonic coupling, such that the 
ground state is always localized for go =/= 0. 



2. Correlation-length exponent 

In the following, we discuss the properties of the 
boundary QPT, controlled by the critical fixed point Cr. 
We start with the correlation-length exponent ^, describ- 
ing the flow away from criticality: The characteristic en- 
ergy scale T* above which quantum-critical behavior is 
observed vanishes as* 



T* cx \t\ v , 



(15) 



where t is a dimensionless measure of the distance to 
criticality, defined such that t > (t < 0) corresponds 
to the localized (delocalized) phase. Upon linearization 
of the RG beta functions around the Cr fixed point, we 
obtain 



4e e 



O (e\e 



2 J2\ 



(16) 



Clearly, v diverges as s — » 1 and r — > 1 together. By 
expanding the square-root in Eq. (|16|) . the inverse cor- 
relation length exponent can be approximated as = 
1 — 2r + s. The same result, valid for small 1 — 2r + s, is 
also obtained in Sec. Ill Cl following an RG expansion valid 
near the strong-coupling fixed point. The divergence of 
v as 1 — 2r + s — > is demonstrated numerically in Sec. 
IIIIB II and the form compared to Eq. (jTHJ) . 
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(a) v 2 A 



Deloc 




( b ) v 2 A 



Deloc 



FIG. 3: Schematic RG flow diagrams for the dissipative 
resonant-level model with p-h symmetry. Although these di- 
agrams are obtained by expansion about the free-impurity 
fixed point and hence are formally valid as r,s — > 1, they 
are confirmed by expansion about the delocalized fixed point 
(Sec. Ill C|l and NRG calculations (Sec. IIII|) to capture the 
correct physics for all r > and < s < I. The hori- 
zontal axis denotes the renormalized bosonic coupling g; the 
vertical axis denotes the renormalized hybridization v. (a) 
r < r+ — (1 + s)/2: Stable fixed points at Deloc and Loc de- 
scribe the delocalized and localized phases, respectively. The 
continuous impurity QPT is controlled by the critical fixed 
point Cr. (b) r+ < r < 1: The delocalized (Deloc) fixed 
point is unstable against finite g. As r — > l - , the Deloc fixed 
point merges with the free-impurity fixed point (FImp). (c) 
r > 1: v is irrelevant. In both (b) and (c), the flow is toward 
Loc for any finite g. 



Near criticality, xi oc (T) is expected to follow a power- 
law form 

for T*«T«T , (18) 

up to a nonuniversal cutoff scale To. This relation defines 
the anomalous exponent r/ x , which governs the anoma- 
lous decay of the impurity "spin-spin" correlation func- 
tion and is calculated via 



Vx = M" 



dlnZ y 



(19) 



v*,g* 



The renormalization factor Z x obeys the exact 
relation* 2 ^ 



z: 



{z g /z f f , 



(20) 



which is graphically represented in Fig.^a). This allows 
us to derive the exact result 



rj x = 2e = 1 — s 



(21) 



at the Cr fixed point, a relation that is borne out by the 
numerical results presented in Sec. IIIII 



4- Conduction electron T-matrix 

The conduction electron T-matrix, describing the scat- 
tering of the c electrons off the impurity, is another im- 
portant observable, being central to the calculation of 
transport properties. For a resonant-level model, the T- 
matrix is given by T(u>) = UqG/(cj) where Gf is the full 
impurity (/-electron) Green's function, graphically rep- 
resented in Fig. QJb). As with the local susceptibility, 
we expect a power-law behavior of the T-matrix spectral 
density near criticality: 



T(uj) oc 



1 

I1-J7T 



for T* < U < T . 



(22) 



3. Response to a local field 

The local impurity susceptibility xioc(T) is the impu- 
rity response to a field applied only to the impurity. 2 
Here, for the spinless resonant-level model under con- 
sideration, the level energy ef plays the role of a lo- 
cal electric field. Defining the impurity "magnetization" 
w; mp = (Sz)i with the pseudospin S z as specified in Eq. 
©, it follows that 



It has been showr>L£ that all critical fixed points for < 
r < 1 in the pseudogap Anderson and Kondo models 
display T[uS) oc |w|~ r as to — * 0, which behavior has been 
observed in a number of separate studies i 33 ' 35 ! 36 

Using the exact relation Zt = Zf/Z%, we can derive 
an exact result for the critical point of the dissipative 
resonant level model: 



Vt = 1 



(23) 



dm-. 



Xloc 



imp 



be j 



(17) 



is nothing other than the impurity capacitance. 



Thus, even though the multiplicative prefactor of the be- 
havior l|2"2"|) is expected to exhibit both r and s depen- 
dence, the power law followed at criticality is identical to 
that of the pseudogap Kondo and Anderson models. 
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(a) 



Xlo 



(b) 



=o 



FIG. 4: (a) Exact relation for the local susceptibility. The 
black triangle denotes the full vertex function and the dashed 
double line denotes the full impurity level propagator, (b) 
The large dot denotes the full hybridization vertex. 



5. Hyperscaling and other critical exponents 

The QCP is expected to satisfy hyperscaling relations 
characteristic of an interacting fixed point, including 
w/T scaling in dynamical quantities. 2 It follows that the 
correlation-length exponent v and the anomalous expo- 
nent r\ x are sufficient to determine all critical exponents 
associated with the application of a local fieldi2iiii For 
example, one can define exponents 7 and 7' through the 
T — > limit of the local susceptibility near criticality: 



X ioc(t<0;T = 0)c<(-i)-T, 7 = 1/(1-^), 
T X XoS >0;T = 0) i = vr hl . 



(24) 



One can also determine critical exponents (3 and 5 asso- 
ciated with the local magnetization mi mp : 



m irop (t>0;T = 0,e / ^0)c<t /3 , /3 = 
m imp (e f ;t = 0,T=0)oc \e f \^ s , S = 2/ Vx -1. 

Thus, near criticality 



(25) 



v/e 2 /4 + 4e(e' - e/2) - e/2 



+ 0(e 2 ,e' 2 ) (26) 



and 
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l + (D(e 2 ,e' 2 ) 



where, in contrast to Eqs. (|2"l"j) and (J33J) , the higher-order 
corrections do not cancel. Section IIIII reports NRG re- 
sults for several of these critical exponents that demon- 
strably obey the hyperscaling relations. 



resulting action isi£ 

S = ^ f(uj n ) [iA sgn(uj n )\uj n /D\ r + iA 1 uj n ] f(w n ) 

r/3 

+ 5 + .g o / dr(//-i)0(r,O), 



(28) 



where the local / fermions are now "dressed" by the con- 
duction lines, 

f 7T7* \ / ITT" \ 

A = TrN v% sec ( y j = T a sec (—) (29) 

is a nonuniversal energy scale, and A± = 1 + O(uq). For 
r < 1, the |w„| r term dominates the / propagator at 
low energies. Then, dimensional analysis of the bosonic 
coupling (here w.r.t. the Deloc fixed point) yields 



[3o] = 



2r - 1 - s 



(30) 



which implies that an RG expansion can be controlled in 
the smallness of 

2i=l-2r + s. (31) 
We introduce a dimensionless coupling according 

/uJ¥nZ n 



9o = A* e Ar 



KnZ 



-9; 



(32) 



and, following the procedure described in Sec. Ill Bl we 
find that the only contribution to Z g is that shown in 
Fig- [He), which reads (note that Zf = 1 to this order) 



1 + esc 



firs\ g 
V 2 / e 



(27) The RG beta function for g is 



(33) 



(34) 



It is clear from Eq. that for s > and e > 0, there 
exists a critical fixed point at 



P(g) =£9-2csc( — ;// 



C. RG expansion around the delocalized fixed 
point 

In addition to the RG expansion for r — > 1 and s — > 1, 
as described in Sec. Ill Bl a second epsilon expansion can 
be performed around the Deloc fixed point. 

1. RG equations 

To begin, we integrate out the conduction electrons, 
which is an exact operation for the present model. The 



»"-!-(=). (») 

which controls the delocalized-to-localized transition. 
The RG flow diagram is sketched in Fig. [5J 

Note that the critical coupling g* approaches zero as 
e — > + and/or as s — > + , suggesting that beyond 
these limiting cases the delocalized fixed point is un- 
stable towards the localized fixed point. The same in- 
stability has already been deduced for e < [i.e., for 
r > r+ = (1 + s)/2], based on expansion about the free- 
impurity fixed point (see Sec. Ill Bp . The behavior for 
s < is analyzed in the next section. 
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FIG. 5: RG flow diagram of the dissipative resonant-level 
model near the delocalized (Deloc) fixed point for s, e > 0. 
The two stable phases are governed by the delocalized (g — 0) 
and localized (g — oo) fixed points, separated by the critical 
fixed point [g — g* specified in Eq. (f35|) ] . 



2. The regime s < 

For s < 0, the perturbation theory described in Sec. 
Ill CI is singular due to the divergent DOS in the bosonic 
propagator. In this range of s, the delocalized fixed point 
is always unstable against any infinitesimal bosonic cou- 
pling go, which favors the localized fixed point. 

We can gain a better understanding of this instability 
by considering the local bosonic propagator G^^iuin) — 
S<j G<j>(Q> i^n) in the presence of the impurity. Including 
impurity effects via the boson self-energy, the local boson 
propagator is given by 



< + sA s - gl for s > 0, 
uj- s - gl for s < 0, 



(36) 



where A is a momentum cutoff energy scale. Let us dis- 
cuss s > first. For sA s > > 0, the local boson 
propagator is massive, meaning that the ground state for 
the bulk is just the empty state. For g^ > sA s > 0, by 
contrast, the local boson propagator has "negative mass" , 
as a consequence of which the local boson condenses at 
zero temperature with an expectation value (</>o) =/= 0. 
This drives the system to the localized phase where the 
pseudospin operator S z also assumes a nonzero expecta- 
tion value. This reasoning supports the existence of a 
QPT for s > 0, with criticality reached at g* 2 — sA s . 
For s < 0, the local boson propagator always has 
a negative mass, i.e., the impurity is localized. (Techni- 
cally, the impurity induces a bound state in G<f, .) The 
observation that the ground state is always localized for 
s < is consistent with previous studies of the spin-boson 
mode l 16 ' 20 and the Bose-Fermi Kondo model ) 28 ' 29 which 
belong to the same universality class as the dissipative 
resonant-level model in the metallic limit r = 0. 



3. Phase diagram 

The RG flow allows us to deduce that the qualitative 
phase diagram of the dissipative resonant-level model in 
the parameter space specified by r and s is as shown in 
Fig. [TJ The solid line denotes the locus of points satis- 
fying 1 — 2r + s = 0. In the unshaded region to the left 
of the line [i.e., for max(0, 2r — I) < s < 1, or equiva- 
lently \ < r + < r < 1 with r + defined in Eq. {H}], the 
RG expansion predicts a continuous QPT as vq and go 



are varied. For s < max(0, 2r — 1) (shaded area), the 
ground state of the model is always localized for any fi- 
nite bosonic coupling go- This is consistent with the RG 
flow diagrams presented in Fig. [31 where the RG expan- 
sion is carried out for r, s — * f . The phase diagram is 
confirmed by NRG results in Sec. IIIII 



4- Critical exponents 

By linearizing the RG equation around the fixed point, 
the correlation-length exponent at the critical point g* is 
found to satisfy 



- = 21 
v 



;2\ 



(37) 



For the anomalous exponent rj x associated with the 
local susceptibility [Eq. ([18])]. we again have the exact 
property Eq. (|20|) [see also Fig. BJa)], from which it fol- 
lows that 



>h 



1 



(38) 



The exponents (3 and 5 can be obtained from the hyper- 
scaling relations (|25| : 



4(= 



Oil 1 



and 



1 



+ 0(e~ 2 ) 



(39) 



(40) 



The exponent rjx, associated with conduction-electron T- 
matrix, is also found to obey t)t = \ — r [see Eq. ([23])]. Of 
course, all critical exponents for the two RG expansions 
(one for r, s — > 1 and one for 1 — 2r + s — > 0) are expected 
to be compatible since the expansions describe the same 
QPT. In the limit r,s — ► I, the square root of Eq. (fl6)) 
may be expanded to yield Eq. ([3T| . The equivalences of 
Eqs. dHJ) and §§§ for (3 and of Eqs. (J2ZJ) and CP for 5 
are also readily verified. 



III. NUMERICAL RENORMALIZATION 
GROUP 

The NRG method^ has recently been extended to pro- 
vide nonperturbative results for the Bose-Fermi Kondo 
mode h 28 ' 29 In the following, we implement the same ap- 
proach for the spinless resonant-level model ([TJ) , which 
also involves both fermionic and bosonic baths. 

There are three essential features of the NRG: (i) The 
energy axis is logarithmically discretized, introducing a 
discretization parameter A. (ii) The Hamiltonian is then 
mapped to a chain form, with the impurity degrees of 
freedom coupled to the first site only of one or more 
tight-binding chains, (iii) Owing to the discretization, 
the tight-binding coefficients decay exponentially with 
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FIG. 6: (Color online) (a) The lowest NRG eigenstates En 
vs even iteration number N for (r, s) = (0.85, 0.9), hybridiza- 
tion strength To = 0.1, and a range of dissipation strengths 
Bo - B , c = 0, ±10~ 3 , ±10" 2 . The flows are typical of those 
for max(0, 2r — 1) < s < 1 with < r < 1. The levels at the 
critical coupling Bo = Bo, c ~ 0.3731 are shown as bold dotted 
lines while those nearby in the delocalized (Bo < Bo, c ) [local- 
ized (Bo > Bo.c)] phase are shown as solid [dashed] lines. As 
Bo approaches Bo, c in either phase, the levels follow those of 
the unstable critical fixed point down to progressively lower 
temperatures, before crossing over to the levels characteristic 
of the delocalized or localized stable fixed point, (b) NRG 
level flows for (r, s) = (0.975,0.9). In this case, and more 
generally for a < max(0, 2r — 1) with < r < 1, the flow is 
towards the localized fixed point for any Bo > 0, but follows 
the delocalized fixed point down to progressively lower tem- 
peratures as Bo is reduced towards zero. The solid lines show 
the flow for Bo = 0. 



increasing chain length. This allows the problem to be 
solved in an iterative fashion, diagonalizing progressively 
longer finite-length chains and thereby including expo- 
nentially smaller energy scales, Tn ~ DAr N / 2 , at each 
iterative step iV = 0, 1, 2, .... The RG transformation 
relating the effective Hamiltonians at consecutive itera- 
tions eventually reaches a scale-invariant fixed point that 
determines the low-temperature properties of the system. 

In all applications of the NRG, the maximum number 
N a of many-body eigenstates retained from iteration N to 
form basis states for iteration N+l must be truncated for 
sufficiently large N due to the limitations of finite com- 
putational power. The presence of one or more bosonic 
chains introduces additional considerations. First, the 
bosonic Hilbert space must be truncated even at itera- 
tion N = 0, allowing a maximum of Nj, bosons per site 
of a bosonic chain. Second, for problems involving both 
fermionic and bosonic chains, the fact that the bosonic 
tight-binding coefficients decay as the square of those for 
fermionic chains must be reflected in the specific iterative 
scheme employed. That is, only (bosonic and fermionic) 
excitations of the same energy scale should be considered 
at the same iterative step. Thus, while the fermionic 
chain is extended at each iteration, the bosonic chain is 
extended only at every second iteration. These issues, to- 



gether with further details of the implementation of the 
Bose-Fermi NRG, are discussed in detail in Ref. H^. 

The NRG method has provided a comprehensive nu- 
merical account of the quantum-critical properties of a 
number of impurity problems, e.g., the fermionic pseudo- 
gap Kondo and Anderson models, the spin-boson model, 
and the Bose-Fermi Kondo model. In all cases it is 
found that the critical properties (such as exponents) are 
insensitive to the discretization parameter A and con- 
verge rapidly with the number of retained states N s . For 
models involving bosonic baths, critical exponents also 
rapidly converge with increasing bosonic truncation pa- 
rameter Nb. In the following we take A = 3, with all data 
suitably converged for the choice N s — 500 and Nb = 8. 
For convenience we set D = ujq = 1. 



A. Phase diagram 

Figure [6] shows the flow of the lowest NRG eigenstates 
Em of the effective Hamiltonian Hn at even iteration 
numbers N for two representative cases for s > 0: (a) 
1 — 2r + s > and (b) 1 — 2r + s < 0. Figure [6][a) shows 
data obtained for (r, s) = (0.85,0.9) and r = 0.1. Here, 
and for any 1 — 2r + s > 0, the flow is schematized by 
Fig. [3ja), which follows from the perturbative analysis. 
For Bo < Bq c , the NRG flow is towards the delocalized 
fixed point, where the spectrum coincides with that for 
coupling Bo = to the bosonic bath. For Bq > Bq c 
the NRG flow is towards the localized fixed point, where 
the spectrum coincides with that for coupling To = to 
the fermionic band. For B close to B , c , as considered 
in Fig. [HJa), the flow in either case is first towards the 
critical spectrum. The departure from the critical flow, 
at a crossover scale T* that vanishes at Bo = £?o,c, is 
governed by the correlation-length exponent discussed in 
Sec. HUB II 

Figure Efb) shows NRG level flows for (r, s) — 
(0.975, 0.9) and T = 0.1. These flows are typical of those 
for any 1 — 2r + s < and correspond to the perturbative 
RG flows of Fig. Elb) . The localized ground state obtains 
for any Bq > 0. As Bo is reduced towards zero, the levels 
follow those of the delocalized fixed point (obtained for 
Bo = 0) down to progressively lower energy scales. 

Figure [7| shows the phase diagram of the model on the 
r-Bo plane for three different combinations of the bosonic 
bath exponent < s < 1 and the hybridization strength 
r . For all s and r pairs considered, the phase-boundary 
value of Bo decreases monotonically with increasing r 
from that found for a metallic conduction band (r = 0). 
This is particularly clear from the data set obtained for 
s = 0.8 and r = 10~ 3 (circles in Fig. [7|), where the 
metallic system undergoes a continuous QPT at a critical 
Bo, c (r = 0) rs 0.699. With increasing r, and hence grow- 
ing depletion of the conduction electron density of states 
around the Fermi level, the critical dissipation strength 
Bo.c required to localize the system is reduced, as ex- 
pected on physical grounds. Bo iC (r) is found to vanish 
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FIG. 7: (Color online) Phase diagram in the r-Bo plane, 
obtained using NRG for the fixed bosonic bath exponent s 
and the hybridization strength To shown in the legend. For 
< r < r+ = (1 + s)/2, we find a continuous QPT between 
delocalized (Deloc) and localized (Loc) phases. The critical 
dissipation strength Bo iC is found to vanish continuously at 
r = r+. For r > r+ only the localized phase can be ac- 
cessed for Bo > 0. The inset shows the vanishing of £>o,c with 
decreasing 1 — 2r + s in each case, compared to the results 
obtained from the perturbative analysis. 



continuously at r = r + , with r + as defined in Eq. (|14[) . 
This vanishing is illustrated in the inset to Fig. [Jj which 
shows -Bo,c vs 1 — 2r + s on a logarithmic scale. 

For r > r + , localized solutions are found for arbitrarily 
small dissipation strength Bq > 0. The symbols at the 
largest r (= r + ) in each case, which lie at Bq = 0, mark 
the point at and above which no delocalized solutions 
can be found with Bq > 0. Thus, we find that we can 
tune the system to a QPT if, and only if, < s < 1 and 
< r < r_|_, in complete agreement with the scenario 
deduced via the perturbative analyses and illustrated in 
Fig.ffl 

For < r < 1 and s = 1 , we find a line of Kosterlitz- 
Thouless-like transitions between delocalized and local- 
ized ground states, and for s > 1 only the delocalized 
phase is accessed (provided r > 0). For r > 1 and 
s > 1, the essential physics is controlled by the free- 
impurity fixed point, regardless of the couplings To and 
Bq. 

For a given (r, s) pair that exhibits a continuous QPT, 
the critical dissipation strength Bq :C varies with the hy- 
bridization strength Tq as 



,(l-a)/(l- 

o 



(41) 



provided that all scales are small compared to the cutoffs. 
This result, which follows from dimensional arguments 
[Eq. (|4"Tj) can readily be obtained using Eq. (fT5)) ] and is 
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FIG. 8: Critical dissipation strength £>o,c vs hybridization 
strength To for the (r, s) pairs specified in the legend. We 
find that B , c oc Tg, with x = (1 - s)/(l - r). 



as the 



confirmed numerically in Fig. [51 identifies Tq^ 1 
tunneling amplitude analogous to A of the spin-boson 
model, where^ the critical dissipation strength is a c oc 
A 1-5 . A similar result for the Bose- Fermi Kondo model 
finds Bq iC oc Tjf s , with Tk the bare Kondo temperature 
serving as a tunneling amplitude between impurity spin 
states,^ 

It is interesting to compare the location of the phase 
boundary obtained using NRG with that inferred from 
analytical expansion. We have in mind fixing the hy- 
bridization strength Tq and the bosonic-bath exponent s 
(as in Fig. [7]), and finding the critical coupling Bq as a 
function of the conduction-band exponent 7'. However, 
an analysis of the expansion around the free-impurity 
fixed point (Sec. Ill Bp reveals no simple analytical ex- 
pression for the phase boundary, due to the fact that 
the problem is described by a two-parameter flow, which 
cannot be linearized in general. We have therefore ana- 
lyzed the coupled differential flow equations numerically. 
The phase boundary can be obtained by determining the 
eigenvalues and eigenvectors of the linearized RG equa- 
tions near the critical point and then following the RG 
flow backwards along the separatrix. 

The inset of Fig. [7] compares phase boundaries deter- 
mined via NRG (symbols) with those obtained via the 
perturbative RG equations (fT5j) (dashed lines). For the 
range of 1 — 2r + s considered by NRG, Bq c appears to 
vanish as a power law, with an exponent that depends on 
both the bosonic bath exponent s and the hybridization 
Tq. This apparent power law does not reflect the asymp- 
totic behavior, revealed by the perturbative calculations 
to be Bo, c oc e as e — > 0. (This regime is inaccessible 
to NRG because the merging of the critical and delocal- 
ized fixed points with decreasing e make it impossible to 
reliably determine the critical coupling Bq jC .) Neverthe- 
less, we find the level of agreement remarkable and stress 
that there is no fitting procedure involved in making this 
comparison. 

From the expansion around the delocalized fixed point 
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FIG. 9: Crossover scale T* vs \t\ = \B — B 0:C \/B 0]C for the 
(r, s) pairs specified in the legend. In the vicinity of the tran- 
sition (\t\ <C 1), T* oc \t\" . The correlation-length exponent 
v(r, s) is independent both of the hybridization To and of the 
phase from which the QCP is approached. 



fSec. MCj) . where we have a one-parameter flow, it seems 
possible to obtain an analytical expression for the phase 
boundary. However, we have to keep in mind that the 
dressed / propagator in Eq. f2"8"|) contains terms with 
different frequency dependencies, and is dominated by 
\uj n \ r in the low-energy limit only. (The coefficient A\ 
is nonzero in general, except right at the Deloc fixed 
point.) The interplay of the \oj n \ r and u) n terms intro- 
duces a nonuniversal crossover scale into the problem, 
and a proper treatment including elevated energies would 
require a multistage RG scheme, which is beyond the 
scope of this paper. 



B. Critical exponents 
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FIG. 10: (Color online) (a) Correlation-length exponent v vs 
conduction-band exponent r for two values of the bosonic- 
bath exponent s, as shown in the legend. The symbols 
show NRG data, while the dashed lines are the correspond- 
ing perturbative results [Eq. ((16}], expanding about the free- 
impurity fixed point. We find that v~ l vanishes at r = r+, 
in keeping with the qualitatively distinct behavior for 2e = 
1 — 2r + s ^ 0. (b) The same data plotted vs 21. For small 
e, v' 1 « 21, a result [Eq. (|37[) ] (shown as a dotted line) 
obtained by a perturbative expansion about the delocalizcd 
fixed point. 



universality class). By increasing r we find that u(r,s) 
diverges as r — ► r + from below, i.e., as 1 — 2r + s — ► + . 
The dashed lines are the corresponding perturbative re- 
sults [Eq. HH])], with which there is excellent agreement 
for r approaching r + . Figure fTUTb') shows the same data 
plotted vs 2e = 1 — 2r + s. With decreasing e > 0, the 
curves approach the result v~ l w 2e (shown as a dotted 
line), as obtained in Sec. II. C. 3 by an expansion about 
the delocalized fixed point. 



Correlation-length exponent 



2. Response to a local field 



The correlation-length exponent v defined in Eq. JT| 
is readily extracted from the crossover scale T* cx A~ N / 2 
in the NRG level flows between the unstable and either of 
the stable fixed points. Here, N* denotes the NRG iter- 
ation number at which crossover is observed in a chosen 
NRG eigenvalue En ■ (See Refs. [H andU^ for further de- 
tails.) Figure^] shows T* vs \t\ = \B - B , c \/Bo, c for the 
(r, s) pairs specfied in the legend. The dashed lines are 
linear fits to the log-log data, which yield the correlation 
length exponent v(r, s), independent of the hybridization 
strength Tq and the phase (Deloc or Loc) from which the 
QCP is accessed. 

The r dependence of the correlation-length exponent is 
demonstrated in Fig. [TUT a) for two values of the bosonic 
bath exponent s. As anticipated, for r = we find that 
within our estimated numerical error of about 1%, ^(0, s) 
is in essentially exact agreement with v(s) for the spin- 
boson mode l 16 i 20 (and the Ising-symmetry Bose-Fermi 
Kondo model, demonstrated in Ref. [2§ to share the same 



As discussed in Sec. Ill B 31 the response to a field ap- 
plied only at the impurity provides a useful probe of the 
locally critical properties of the model. The inset to Fig. 
[TTTa) shows m imp (t;T — 0) vs t = (B a — B , c )/B _ c for 
(r, s) — (0.85,0.9) and hybridization strength T = 0.1. 
Behaving as a suitable order parameter for the problem, 
Mm P (t;T — 0) is finite in the localized phase (i > 0), 



0) 
0^ 



for t 3> 1 and van- 
In the delocalized phase 
0) = 0. The main part of Fig. 
0) vs t > on a logarithmic scale, 



saturating to mi mp (i;T 
ishing continuously as t 
(t < 0), {t-T 
[TT] shows m; m p(t;T 
from which the power-law behavior Eq. (|25p is clearly ap- 
parent. The exponent (3 is found to be (3 — 0.601(2). At 
the QCP (t = 0), the dependence of m imp (t = 0,T = 0) 
on the field ej defines the exponent S according to Eq. 
(f2"S")) . We typically observe such power-law behavior over 
several orders of magnitude of £/, as shown in Fig.fTlTb'). 
For (r,s) = (0.85,0.9), 1/5 = 0.052(1). 

We note that for < r < 1 and s = 1, mi mp (i;T = 
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FIG. 11: (Color online) Critical exponents /3 and S, defined 
in Eq. J25}, for (r, s) = (0.85, 0.9) and To = 0.1, where B , c « 
0.3731. (a) Continuous vanishing of order parameter mj mp vs 
t = (Bo-~Bo >c )/Bo, c as t — > + with characteristic exponent /3 
(extracted as the limiting slope of the data on a logarithmic 
scale). The inset shows the data on an absolute scale, (b) 
Variation of mi mp (T = 0) with local level energy Ef at the 
critical point t = 0. The data clearly follow a power law for 
small Ef, defining the exponent 8. 
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FIG. 12: Static local susceptibilty Xl oc (T) vs T for (r, s) = 
(0.2,0.5), To = 0.1, and B = 0.4902 (circles), 0.5002 as fl ,c 
(stars), and 0.5102 (squares). The anomalous exponent in the 
quantum-critical regime is found to be r\ x = 1 — s, independent 
of r. See text for further discussion. 



0, Sf — + ) undergoes a jump at the critical point t = 0. 
Here, the essential behavior has been discussed in Refs. 
SHallS and [13 for the case (r, s) = (0, 1) relevant to 
charge fluctuations on a metallic island subject to elec- 
tromagnetic noise. 

We calculate the static local susceptibility via 



Xio C (T) 



lim 



in 



imp 



(42) 



(T = 0) van- 



In the delocalized phase B < B c 
ishes linearly with Ef and thus Xioc(T) w const, for 
T <C T* . In the localized phase Bq > Bq jC , m; mp is 



nonzero as e f -> with xi oc (T < 7 1 *) cx 1/T. In the 
quantum-critical regime T* -c T -C T , XiocC^) diverges 
as a power law with an anomalous exponent 7y x defined 
in Eq. (fT5|) . For all (r, s) pairs considered (such that 
1 — 2r + s > and a critical fixed point exists), we find 
that 



% = 1 - s - 



(43) 



independent of r. The behavior described above is clearly 
illustrated in Fig. IT2J which shows three data sets for 
(r,s) — (0.2,0.5): one at the critical coupling and one 
close to it in either phase. In this example, we extract 
r] x = 0.499(2). 



3. Hyperscaling 



As discussed in Sec. Ill B 51 critical exponents for the 
present model are expected to obey hyperscaling rela- 
tions derived via a scaling ansatz for the critical part 
of the free energy that assumes the critical fixed point 
is interacting^ This expectation is borne out by the 
numerical analysis: we find hyperscaling relations to 
be obeyed to within the estimated error (typically less 
than 1%) across the range of (r, s) displaying critical be- 
havior. For example, for the case (r, s) — (0.85,0.9), 
\jv = 0.082(1) and J] x = 0.101(2). Thus, the values 
/3 = 0.601(2) and 1/6= 0.052(1) extracted from the data 
presented in Fig. [TT]obey Eqs. (f2"S"|) to within numerical 
uncertainty. 



C. Spectral function 

We now turn to the single-particle spectral function 
A(uj), calculated via 



m )\ 2 ~ 7p - E n + E m ), 

(44) 

where \m) is a many-body eigenstate of NRG iteration 
N, and Z — J2 n ex P( — PE n ) is the partition function; 
A(u>) = A(— oj) for the p-h symmetric parameters stud- 
ied. The discrete delta-functions are Gaussian broadened 
on a logarithmic scale: a standard NRG procedure dis- 
cussed, e.g., in Ref. IH. We set the broadening parame- 
ter b such that A(uo) for the simplest resonant-level model 
(with r — 0, Bq — 0, and ej — 0) is in optimal agreement 
with the exact result A(w) = 7r _1 r /(a; 2 + Tq). 

Figure fT3T a) shows A(ui) vs \u>\ on a logarithmic scale 
for r = 0.65, s — 0.8, To = 10~ 3 , and the dissipation 
strengths B < B 0tC specified in the figure caption. For 
the delocalized phase Bq < £?o,o we find that the dissi- 
pation does not alter the asymptotic low-frequency be- 
havior of A(u>) found for B = 0, i.e., 
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FIG. 13: (Color online) (a) Spectral function A(lj) vs for 
r — 0.65, s = 0.8, To = 10 -3 , and three values of the dis- 
sipation strength: Bo = (dotted line), Bo = Bo, c — 10~ 5 
(solid line), and B = B 0>c = 0.03247113 (thick dashed line). 
At the quantum critical point Bo = Bo, c , A{u) oc \u\~ r , 
which behavior is also followed for Bo close to Bo, c and 
[w| S> T*. In the delocalized phase Bo < Bo, c , there is a 
crossover in A(ui) to the behavior Eq. (|45p for \u\ <C T„. 
For the data shown, T* ~ O(10" 26 ). (b) The crossover 
behavior is more readily seen in the modified spectral func- 
tion J-(lo) = 7rFo sec 2 (-|r) \ui\ r A(uj), which shows the ultimate 
low-o; behavior T{lo = 0) = 1 throughout the delocalized 
phase Bo < Bo, c , < T(w = 0) < 1 for Bo = Bo, c , and 
T(oj — 0) = throughout the localized phase Bo > Bo, c - In 
order of decreasing crossover scale, delocalized-phase spectra 
are shown for Bo = (dotted line) and for Bo, c — Bo = 10 -3 , 
10 -4 , 10 -5 , and 10 -6 (solid lines); localized-phase spectra 
(dashed lines) are shown for Bo = 0.05 and Bo — Bo. c = 10~ 3 , 
10 -4 , 10 -5 , and 10~ 6 . The critical spectrum is shown as a 
thick dashed line. 

For Bq = the spectrum is identical to that obtained for 
the noninteracting (U = 0) limit of the (spinful) pseudo- 
gap Anderson model at p-h symmetry, where the result 
Eq. (|35]l holds for < r < lJ^ Moreover, it is known2iji2£ 
that the form Eq. (|¥5|) persists throughout the Kondo- 
screened phase of the pseudogap Anderson model with 
interactions present (i.e., for all U < U c ), which in the 
p-h symmetric case is confined to < r < i. 
In the vicinity of the QCP, Bq w i?o,c we find 

A(u) = cos2 (y) M~ r for T * « 1^1 « T °' 

(46) 

where c(r, s) < 1 and To is a high-frequency cutoff set 
by the bare hybridization strength Tq. This behavior 
confirms Eqs. ((22) and (|2lij). 

In the localized phase, by contrast, A{uS) vanishes as 



lj -> 0: 

A(w)oc|w| a for|w|<T*. (47) 

The exponent a is positive, and in general depends on 
both r and s. 

The crossover between these behaviors is more read- 
ily apparent in the modified spectral function T{lo) = 
7tFo sec 2 (7rr/2)|u;|'"A(a;). Any low-frequency divergence 
of A(uj) is canceled in T{tjj), and .F(O) = 1 is pinned 
throughout the delocalized phase of the model. As 
discussed in the context of the pseudogap Anderson 
model ; 34 ' 35 ! 36 this generalizes the well-known pinning 
iiTqA({)) = 1 of the spectral function for the regular 
(r = 0, fermionic) Anderson model. In the delocalized 
phase, the scale T*, playing the role of a renormalized 
tunneling amplitude, is then manifest as the width of the 
pinned resonance at the Fermi level lo = 0, vanishing as 
B ° B o,c- 

Figure fT3fb) shows !F{ui) vs for r = 0.65, s — 0.8, 
r = 10~ 3 , and the B values specified in the figure cap- 
tion. Throughout the delocalized phase (0 < Bq < Bq iC ), 
^-"(0) = 1 remains satisfied to within a few percent, as 
is typical for NRG. Close to the QCP in either phase, 
!F{ui) w c(r, s) down to the scale T*. 

We close by considering the single-particle spectrum 
for the case of a metallic fermionic density of states 
(r = 0) and Ohmic dissipation (s = 1). Here the model 
describes charge fluctuations on a quantum dot or res- 
onant tunneling device close to a degeneracy point and 
subject to electromagnetic noise. The essential physics — 
a Kosterlitz-Thouless-like QPT between delocalized and 
localized states — has been investigated in a number of 
earlier studie o 29 ' 30 i 37 ' 38 i 39 , e.g., via a Bose- Fermi Kondo 
model, and we will not repeat the discussion here. We 
simply show, in Fig. 1141 the spectrum for Tq — 0.001 and 
a range of dissipation strengths; for Bq = 0, A(ui) is of 
Lorentzian form. The vanishing width of the central res- 
onance as Bq — > Bq c indicates a suppression of tunneling 
between dot and leads due to the noisy electromagnetic 
environment. 



IV. CONCLUSIONS 

In this paper, we have analyzed the phase diagram and 
the quantum phase transitions of a paradigmatic quan- 
tum impurity model with both fermionic and bosonic 
baths, namely a dissipative resonant-level model. For 
weak dissipation, the resonant tunneling of electrons is 
renormalized due to the friction of the bosonic bath, but 
the ground state remains delocalized. For strong dissi- 
pation, by contrast, the tunneling amplitude renormal- 
izes to zero in the low-energy limit leading to a local- 
ized ground state. We have employed both analytical 
and numerical techniques, utilizing epsilon expansions re- 
cently developed in the context of the pseudogap Ander- 
son and Kondo model, and an extension of Wilson's nu- 
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FIG. 14: (Color online) ttToA(lu) vs u) for the case of a metal- 
lic fermionic density of states (r = 0) and Ohmic dissipation 
(s = 1) for To = 0.001. Spectra are shown for increasing dissi- 
pation strength Bo (see legend) in the delocalized phase. The 
spectrum is a simple Lorentzian for Bo — 0, and the vanish- 
ing width as Bo — * Bq c indicates a suppression of tunneling 
between the local level and the conduction band. 



presented numerical results, making contact with earlier 
bosonization studies of related models. 

We finally mention a few applications. In the context 
of nanostructures, a resonant-level model may describe 
the tunneling of electrons between a lead and a small 
island or quantum dot. 40 i 42 Taking into account elec- 
tromagnetic noise of a fluctuating environment directly 
leads to a model of type Jl}, provided that the spin de- 
gree of freedom of the electrons can be neglected (e.g., if 
electrons are spin-polarized due to a large applied mag- 
netic field). Related situations, mainly corresponding to 
bath exponents r = and s = 1, have been discussed in 
the literature i 30 i 39 Apart from the common situation of 
ohmic noise (s = 1), sub-ohmic dissipation (s < 1) can 
occur, e.g., in RLC transmission lines which display a y'u; 
spectrum in the R-dominant limits Further, a bath with 
r = 1 may be realized using Dirac electrons of graphene 
or quasiparticles of a d-wave superconductor. 
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